Novel anisotropy in the superconducting gap structure of Bi 2 Sr 2 CaCu208+5 probed by 

quasiparticle heat transport 
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Since the nature of pairing interactions is manifested in the superconducting gap symmetry, 
the exact gap structure, particularly any deviation from the simple d x i_ y 2 symmetry, would help 
elucidating the pairing mechanism in high-T c cuprates. Anisotropic heat transport measurement in 
Bi2Sr2CaCu208+i reveals that the quasiparticle populations are different for the two nodal directions 
and thus the gap structure must be uniquely anisotropic, suggesting that pairing is governed by 
interactions with a rather complicated anisotropy. Intriguingly, it is found that the "plateau" in the 
magnetic-field dependence of the thermal conductivity is observed only in the 6-axis transport. 

PACS numbers: 74.25.Fy, 74.20.Rp, 74.72.Hs 



It has been established Q that in the high-T c cuprates 
the superconducting gap A has essentially the d x 2_ y i 
symmetry, A(k) = Ao cos 20, which is depicted in Fig. 1 
(Ao is the maximum gap and 4> is the azimuthal angle). 
Although this information is very important, it does not 
help much in clarifying the high-T c mechanism, because 
all of the proposed theories of the high-T c superconduc- 
tivity have successfully adopted the d-wave symmetry. It 
is thus important to obtain more detailed information on 
the gap structure, because any deviation from the simple 
d x 2_ y 2 symmetry would point to the role of a particular 
interaction in the paring. However, not much is known 
about the precise gap structure of the high-T c cuprates 
PL partly because the experimental techniques to selec- 
tively probe the nodal structures are limited ||. 

The most prominent feature of the d-wave gap is that 
it has four nodes where the gap magnitude vanishes. At 
finite temperature below T Cl quasiparticles (QPs) are al- 
ways thermally excited near these nodes from the super- 
conducting condensate; how many QPs are excited de- 
pends on the "slope" of the gap at the node j| , because 
a steeper gap allows less phase space available below a 
given energy of ksT. Since these QPs carry heat (in con- 
trast to the Cooper pairs that do not carry heat) , one can 
investigate the gap structure of unconventional supercon- 
ductors by measuring the QP heat transport in the su- 
perconducting state p. Moreover, because the QPs can 
only travel along the nodal directions Q, one can gain 
information about different nodes by tuning the direction 
of the heat current to the respective nodal directions; this 
makes it possible to study the detailed anisotropy of the 
gap structure in high-Tc cuprates. 

In recent years, the heat transport has been conve- 
niently used to probe the gap structure of the high-T c 
cuprates. For example, it was reported || M that the in- 
plane thermal conductivity showed a four-fold-symmetric 
change when the heat current was fixed along an antin- 
odal direction and an applied magnetic field was 360° ro- 
tated in the Cu02 plane, giving maxima for fields parallel 
to the nodal directions. This feature helped establishing 
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FIG. 1: Schematic pictures of the d x 2_ y 2 -wave gap and the 
present experiment, (a) Shape of the gap, A(k) = Aocos20, 
is plotted together with the directions of the a and b crystallo- 
graphic axes of BSCCO. (b) Structure of the two-dimensional 
CuC>2 plane, (c) Two samples used for the measurements of 
the a- and 6-axis transports are essentially the same crystal. 



the d x 2_ y 2 symmetry of the superconducting gap. The 
heat transport measurement was also employed in an in- 
triguing proposal of a field-induced transition in the gap 
structure; namely, Krishana et al. reported that at 
low temperatures the magnetic-field (if) dependence of 
the thermal conductivity of P^Si^CaC^Os+a (BSCCO) 
showed a "plateau" in high fields, and they proposed that 
this plateau region may correspond to a fully-gapped 
state where there are no QPs to carry heat. This ex- 
periment attracted a lot of attention because the field- 
induced phase transition in the superconducting state has 
a significant implication on the mechanism of the super- 
conductivity; however, the origin of the plateau, as well 
as the reproducibility of the experimental result, are cur- 
rently under intense debate (see Ref. || and references 
therein). In this Letter, we report a striking in-plane 
anisotropy of the QP heat transport, which not only sorts 
out the "plateau" issue but also suggests an unexpected 
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FIG. 2: Anisotropic transport in zero field, (a) Resistiv- 
ity shows a very small, if any, in-plane anisotropy, indicating 
that the electronic properties are essentially isotropic in the 
normal state, (b) Thermal conductivity k shows a peculiar 
in-plane anisotropy below T c that suggests a difference in the 
QP populations along the two nodal directions. 



anisotropy in the superconducting gap symmetry. 

We measure the thermal conductivity k of high-quality 
BSCCO single crystals using a conventional steady-state 
technique |5| ^] . The crystals used here are grown by the 
floating-zone (FZ) technique [|ll| and are annealed to be 
slightly underdoped (zero-resistance T c is 90.5 K); their 
crystallographic axes are determined by the X-ray Laue 
analysis within the accuracy of a few degrees. To obtain 
reliable data on the in-plane anisotropy, we cut a piece 
of high-quality single-domain crystal into three parts and 
measure the a and b axis transports separately [Fig. 1(c)]; 
both the resistivity (p a and pb) and the thermal conduc- 
tivity (n a and Kb) are measured on the same samples. 
Note that in BSCCO the a and b crystallographic axes 
are along the nodal directions [Fig. 1(a)]. The mag- 
netic field is applied perpendicular to the ab plane and 
the independence of n is measured with the field-cooled 
procedure |^] to avoid complications that are associated 
with the vortex-pinning-related hysteresis fL2[ . 

Figure 2 displays the temperature dependences of p 
and k for the a and b directions; the resistivity shows 
only a small, if any, in-plane anisotropy, while the ther- 
mal conductivity demonstrates a non-trivial anisotropy 
that cannot be attributed to an error in the geometrical 
factors; namely, the peak height at ~70 K is smaller along 
the a axis and yet K a becomes larger than Kb below ^30 
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FIG. 3: Magnetic-field dependences of the anisotropic ther- 
mal conductivity at (a) 12.5 K and (b) 10 K. The "plateau" 
is observed only for the b direction. The data are normalized 
by the values in zero field. (At 12.5 K, n a and Kb are 2.07 and 
1.96 W/Km; at 10 K, «„ and K b are 1.96 and 1.81 W/Km). 



K. This anisotropy in k(T) suggests that the number of 
thermally-excited QPs is larger for the nodes along the a 
axis, because the QP population near the nodes is more 

(The 
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directly reflected in k at lower temperatures 
anisotropy in k above T c is most likely due to phonons; it 
has been estimated that phonons are responsible for 
roughly 2/3 of the total k at T c .) We note that the un- 
certainties in the absolute values in Fig. 2 are less than 
10% for both p and n. 

The magnetic-field dependence of k reveals more strik- 
ing in-plane anisotropy of the heat transport. Figure 3 
shows the result of our precise ii-dependence measure- 
ment (where the relative accuracy is <~0.1% ||) for K a 
and ki, at 12.5 and 10 K up to 6 T. Interestingly, only 
Kb shows the plateau-like feature, indicating that there 
cannot be a true "fully-gapped" state where all the nodes 
are lifted. (If a "full gap" causes the QPs to disappear, 
there should be no QP heat transport in any direction.) 

To further investigate the magnetic-field dependence 
of k down to lower temperatures and up to higher mag- 
netic fields, we have measured K a and Kb in a 3 He refrig- 
erator equipped with a 16 T magnet (with a bit worse 
signal-to- noise ratio). Figures 4(a) and 4(b) show the 
if-dependences of K a and Kb, respectively, down to 0.36 
K and up to 16 T; the trend is clear that K a tends to 
decrease with H while Kb tends to flatten or increase 
with H . Note that, since the phonons contribute to the 
heat transport in parallel to the QPs, the observed ther- 
mal conductivity is a sum of the QP term kqp and the 
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FIG. 4: Magnetic-field dependences of the anisotropic ther- 
mal conductivity at lower temperatures and in higher fields. 
(a,b) H-dependences of K a and Kb at various temperatures 
down to 0.36 K. (c-e) Comparisons of K a and Kb in absolute 
magnitude at representative temperatures, (f) Possible shape 
of the gap suggested by the data; the slope of the gap at the 
fe-axis nodes is steeper than that at the a-axis nodes. 



phonon term n p h- For BSCCO, it has been discussed 
H H |l2) that more than 90% of the total n is due to 
K p h around 10 K and this n p h is independent of H ; thus, 
although the independence of the total k is rather weak 
in Figs. 4(a) and 4(b), kqp is actually changing a lot 
with H. 

The physical origin of this anisotropic behavior can 
be inferred by plotting n a and Kb at the same tempera- 
ture together, taking the absolute values for the vertical 
axis [Figs. 4(c)-4(e)]. In the simple kinetic theory, the 
QP thermal conductivity kqp is proportional to both the 
number of QPs available for the heat transport, ./Vqp, and 
their mean free path Zqp; namely, kqp oc -/Vqp Zqp. In d- 
wave superconductors, the quantized vortices produced 
by applied magnetic fields are expected to both induce 
QPs near the nodes and scatter them |I3, |l]| [l6]. If 
there are already a large number of QPs along the a-axis 
in zero field, the increase in -/Vqp due to the vortices may 
be inconsequential, while the vortex scattering will cause 



a noticeable reduction in Zqp; as a result, kqp along the 
a-axis would tend to decrease with increasing H . On 
the other hand, if there are only a small number of QPs 
along the 6-axis in zero field, both ./Vqp and Zqp along 
the 6-axis will be noticeably changed due to the vortices 
and the changes in Nqp and Zqp will tend to cancel each 
other, leading to a weaker independence along the 6- 
axis. (Actually, one of the probable explanations of the 
"plateau" in k(H) is that it comes from an exact can- 
cellation of the two effects |pj[.]) This phenomenology 
explains the behaviors of K a (H) and Kb(H) in Figs. 4(c) 
and 4(d) (7.5 and 3.0 K) and is consistent with their 
relative magnitude. At lower temperature [Fig. 4(e)], 
it appears that the changes in ./Vqp dominate over the 
changes in Zqp, causing k to increase with H. (Unfortu- 
nately, more quantitative analysis of the data is hindered 
by the lack of proper understanding of the vortex scat- 
tering of the QPs in the cuprates, which remains rather 
controversial |l5|, |l6||.) 

Thus, the novel anisotropics in both the temperature 
dependence and the magnetic-field dependence of the 
thermal conductivity consistently tell us that there are 
more QPs near the nodes along the a-axis; this gives 
strong evidence that the gap structure notably devi- 
ates from the d x 2_ y 2 symmetry. Mixing of a small s- 
wave component, which the orthorhombic distortion may 
cause Q, would not produce the observed anisotropy. 
One simple possibility to make the gap consistent with 
the observed anisotropy is to introduce "higher harmon- 
ics"; for example, for a gap with an expression A(k) = 
A o cos20 — Aisin40, the slope at the nodes is gentle 
along the a axis and steeper along the b axis [Fig. 4(f)] 
[Q. Of course, some non-analytic modulation of the 
d x 2_ y 2 symmetry may also be possible. 

Since the reproducibility of the plateau in k{H) has 
been a matter of debate [||, we comment on the repro- 
ducibility issue. We reported previously that the plateau 
was not observed in the majority of samples 0; recently, 
we found that this majority of samples were cut along 
the a-axis, because our crystals usually grow along the 
a-axis in the FZ furnace. In fact, the k(H) data pub- 
lished in Ref. |9| agree with the n a data reported here. 
To check whether our salient observation reported here 
is reproducible, we have measured another 6-axis sample 
(Sample B2) that was taken from a different batch where 
the crystal grew along the 6-axis. We have also measured 
a sample (Sample M) which was cut along the direction 
45° rotated from the a-axis (Cu-O-Cu bond direction). 
As shown in Fig. 5, Sample B2 (which is optimally-doped 
with T c — 95 K) reproduces the behavior of the first 6- 
axis sample [Fig. 4(b)], and Sample M shows the behav- 
ior which is essentially an average of K a (H) and Kb{H). 
These results give confidence in the reproducibility as 
well as the intrinsic nature of our observation. 

It is useful to note that the angle-resolved photoemis- 
sion experiments have reported almost completely four- 
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FIG. 5: H dependences of « in (a) the second 6-axis sample 
(Sample B2) and (b) a sample cut along the Cu-O-Cu bond 
direction (Sample M). 



fold-symmetric intrinsic Fermi surface for BSCCO [fL9|, 
and therefore the novel anisotropy of the gap found here 
is not likely to be due to a Fermi surface anisotropy 
[p0|| . (In fact, the normal-state resistivity shows negligi- 
bly small in- plane anisotropy [Fig. 2(a)], which is consis- 
tent with the almost-isotropic Fermi surface.) Therefore, 
the striking anisotropy in the gap structure suggests that 
the driving force of the gapping, or the "glue" to bind the 
electrons into Cooper pairs, is anisotropic itself. From 
this point of view, it would be difficult to imagine that the 
antiferromagnetic interaction on the nearly-square Cu02 
planes is wholly responsible for the pairing, though it is 
often considered to be the most likely candidate. On the 
other hand, the phonons are expected to be anisotropic 
in BSCCO because of the pronounced modulation struc- 
ture pi] ) and thus the role played by the phonons in 
pairing is an interesting possibility. Another possibility 
is a role of the charge inhomogeneity [ p2| , which could 
be related to the isotropic or nematic phases of the elec- 
tronic liquid crystal |2^, ^4|, though the impact of such 
inhomogeneous structure on the gapping and on the QP 
behavior is not understood yet. In any case, the novel 
anisotropy in the gap structure suggested by the present 
experiment should help clarifying how the superconduct- 
ing gap is formed in high-T c superconductors. 

In summary, we have measured the thermal conduc- 
tivity k of BSCCO along the a and b axes in essen- 
tially the same high-quality single crystal down to 0.36 
K. It is found that the temperature dependence and the 
magnetic-field dependence of n both show peculiar in- 
plane anisotropy, and the observed anisotropy is consis- 
tent with the interpretation that there are more QPs near 
the nodes along the a axis. This result gives evidence that 
the d-wave gap is uniquely distorted in BSCCO, which 
suggests that the electron pairing is contributed by the 
interaction that develops a pronounced anisotropy be- 
tween the a and h axes. Tn addition, our measurement 



field dependence; this new information should help elu- 
cidating the true origin of this intriguing plateau. 
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